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Abstract 

This work builds on the Volterra series formalism presented in [D. W. 
Dreisigmeyer and P. M. Young, J. Phys. A 36, 8297, (2003)] to model 
nonconservative systems. Here we treat Lagrangians and actions as 'time 
dependent' Volterra series. We present a new family of kernels to be used 
in these Volterra series that allow us to derive a single retarded equation of 
motion using a variational principle. 
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1 Introduction 



The central question addressed in this paper is: 'How can one have a single 
retarded equation of motion arise from the use of a variational principle?' Having 
a dissipative equation of motion arise from the use of a variational principle has a 
rather long history in mechanics. One of the central results in this area is Bauer's 
1931 corollary 0: 

Corollary 1.1 The equations of motion of a dissipative linear dynamical system 
with constant coefficients are not given by a variational principle. 

Bauer's proof of this corollary relies on two implicit assumptions. First, there 
are no extra equations of motion that arise. Also, only integer ordered derivative 
operators are used in the action and Lagrangian. It is the latter assumption that we 
will eventually take advantage of to avoid Bauer's corollary. 

There has been some effort devoted to bypassing Bauer's result. One of the 
earliest attempts was by Bateman ||2|. What he did was allow a dual equation of 
motion to arise. Let us illustrate this via the dissipative harmonic oscillator. If we 
start with the Lagrangian 

C 

L = mxy + —{xy — xy) — mu'^xy, (1) 
where C is a constant, we would have the following equations of motion 

mx + Cx + muj^x = (2) 
my — Cy + mui^y = . (3) 

Equation © is what we want for our harmonic oscillator (i.e., it is retarded or 
causal). Equation © is a time reversed version of ©. Bateman's method is not 
particularly general. Also, the appearance of an advanced (time reversed or anti- 
causal) equation can be considered a drawback to this procedure. Our universe 
appears to be causal so there does not seem to be any compelling reason why 
anti-causal effects should arise from any correct action. 

A rather novel attempt to get around Bauer's corollary was explored by Riewe 
in a series of papers \9 . 10 1 . Riewe tried using fractional derivatives in the action to 
have nonconservative equations arise. In particular, Riewe did not use Bateman's 
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method of having an anti-causal dual equation. However, Riewe's equations of 
motion are acausal. In order to circumvent this, it was suggested that all anti-causal 
operators be replaced with their causal counterparts. (Riewe's method is rather 
involved so we do not provide any illustrative examples.) This procedure does 
allow for much more general equations than Bateman's method. However, it must 
be remembered that the actual equations derived are acausal and the procedure to 
change these into causal equations is rather ad hoc. 

In examining Riewe's approach, Dreisigmeyer and Young showed that his 
procedure of replacing anti-causal with causal operations may not be a wise idea 
[6 1. While still employing fractional derivatives, Dreisigmeyer and Young instead 
allowed for a causal equation and an anti-causal dual equation to arise. Thus, their 
procedure can be considered a generalization of Bateman's. The appearance of 
an anti-causal equation is still very troublesome. Since both 10 [TD1| and ^ deal 
with fractional derivatives, we will later refer to these techniques as fractional 
mechanics. 

There are some other approaches to our basic problem. Tonti fOl and Arthurs 
and Jones [ 1 1 did develop a procedure based on the convolution product that is 
useful for the harmonic oscillator. However, this does not seem to generalize to 
higher ordered potentials. So we need to question if this is the correct procedure to 
follow. Caldeira and Leggett modelled nonconservative systems by coupling 
them to an environment. The environment is modelled as a collection of harmonic 
oscillators which results in the Lagrangian 

n=l J 

where q is the system's coordinate and the g„ are the environment's coordinates. 
Here the combined system and environment is conservative while the system 
alone may be dissipative. This procedure allows for the introduction of very 
general dissipation terms into the system's equation of motion. For example, 
fractional derivatives can be included. However, the microscopic modelling of the 
environment makes (HI) rather complicated. In particular, we would generally have 
a countably infinite system of equations resulting from (21). Also, the question of 
causality is 'put off to the end. That is, we enforce causality in g's equation of 
motion by only using the causal solutions of the g„'s equations of motion. In other 
words, we assume the environment acts causally only when we are solving it's 
equations of motion. While this is not particularly disturbing, it does mean that 
causality is not necessarily enforced in the action. 
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In fF| the idea of treating the actions and Lagrangians as Volterra series was 
introduced. (Volterra series are the generalization to functionals of the power series 
concept.) This is a rather powerful framework. It allows us to build up Lagrangians 
(and hence actions) in a rather systematic way by tailoring the kernels in the Volterra 
series to our requirements. As originally stated, we desire to derive a single causal 
equation of motion for a system. We will accomplish this by choosing the correct 
family of kernels for our Volterra series. Our paper is organized as follows. We 
review the distributional approach to fractional differentiation in section El In 
section El we introduce Volterra series. Our variational principle is developed in 
section HJ A discussion of our results and possible future research follows in 
section |5l 

2 Fractional derivatives 

We will now briefly review fractional derivatives using the distributional approach. 
(A fuller discussion of this material can be found in OEHHIl.) First, define the 
generalized functions 



where r(a) is the gamma function. These distributions will allow us to define 
two different types of fractional derivatives. One of these will be causal while the 
other will be anti-causal. The anti-causal derivatives need to be avoided in our 
equations of motion. 

Left fractional derivatives (LFDs) of order a of a function q(t) are defined by 




(5) 



and 




(6) 



a 



(V) 





(8) 
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where V is the generalized derivative. The LFDs are causal operations. Hence, 
we will only want LFDs to be in our equations of motion. 

Right fractional derivatives (RFDs) of order a are given by 

fitli] ■= *:<,(*)* 9(0 (9) 

where now q{t) = for t > 6. Notice that instead of (EJ), we have 

.DHg] = i-irV-q. (10) 

The RFDs are anti-causal operations. So, we do not want the RFDs to appear in 
our equations of motion. 

The basic problem for fractional mechanics is to avoid the appearance of 
RFDs in the equations of motion. Neither Riewe |l9l [T0| nor Dreisigmeyer and 
Young |6| were able to remove RFDs completely. Riewe suggested replacing 
RFDs with LFDs in the resulting equations of motion. Dreisigmeyer and Young 
suggested letting the LFDs appear in one equation while the RFDs appeared in a 
dual equation. Neither of these methods is entirely satisfactory. We will shortly 
present a formalism that completely avoids RFDs. Before that, let us look at 
Volterra series. 



3 Volterra Series 

Now we are going to develop some of the Volterra series concepts that we will need. 
Our treatment of Volterra series here is somewhat different than that presented in 
|f6|. Here we treat Volterra series as expansions of (generalized) functions that 
depend on time. This is the viewpoint adopted in the nonlinear systems theory that 
is somewhat popular in electrical engineering (see, e.g., ifTDl ^ In the Volterra 
series were treated as expansions of functionals only (i.e., no time dependence). 
This change in viewpoint will make all the difference in achieving our goal of a 
single causal equation of motion. 

For completeness let us first review the Volterra series treatment in fE\, see 
also |12|. (Please note, we will not review the fractional mechanics formalism 
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developed by Dreisigmeyer and Young in this paper since this will be largely 
irrelevant for our current work and also take us too far afield.) In fSl only 
functionals were treated as (time independent) Volterra series. For some functional 
V[q] define the kernels 

Ki^)ir, r) ■= (11) 

Notice that the Kn\-ys are symmetric under an interchange of the Xj's. For 
example, K2 \ti,T2) = K2 \t2, Ti). Now introduce the convenient notation 

i^W^g" := / Kl'\n,...,TMrn)---q{n)dTr,---dn. (12) 

Then we can represent the functional V[q] as the Volterra series 

V[q] = (13) 

n=l 

The V[g] were taken as the actions in f6l. (This is why we can take Kq*^ = V[0] = 
in (fT3t without any loss of generality). This use of Volterra series proves to be too 
restrictive [5 1. So let us use a more 'dynamic' form of the Volterra series concept. 

Instead of expanding the actions in Volterra series like ([T3t . we will now 
expand the Lagrangians in Volterra series like 1.1 IJ 

OO ^ 

^[^J^] = y,-,Kn{T)kq^ (14) 

n=l 

where the ★ notation is as in dT^ and the Knir) will now be of the form 

Knir) = Kn{T,Ti, . . . ,Tn). (15) 

(We assume that Kq{t) = in (fT?b.) We call akemel K„(t, ti, . . . , t„) symmetric 
if it is symmetric under any interchange of the Xj, i = 1, . . . , n. A kernel is called 
stationary if there exists a kernel k„ (^i, . . . , (^„) such that 

Knir-Ti,...,T-Tn) = J'^^ (t, Ti , . . . , T^) (16) 

for all T, ri,...,r„. 
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The r] will generally be r-dependent distributions that also depend on the 
function q. So, in a sense, the Lagrangians can also be viewed as r-dependent 
functional. Hence the notation C[q; t] used in (fT?b . This is also why we allow 
the extra r to appear in (fTSl versus (fTTT) . 

From (fT4b we will form the actions 



The actions given by (flTl are t-dependent Volterra series. So they can be thought 
of as distributions or t-dependent functionals as the Lagrangians were. It is also 
useful to view the actions as the anti-derivatives of the Lagrangians. 

Our entire problem can now be succinctly stated: 'Find the correct kernels 
Kn{T, Ti, ... ,Tn) to usc in (fT?b .' The Volterra series framework helps clarify the 
problem of deriving a single causal equation of motion. It reduces the whole 
question to that of finding a correct family of kernels for our Volterra series. Our 
attention now turns to this problem. 

4 Causal fractional Lagrangian mechanics 

Systems and control theorists tend to be obsessed with the question of causality. So 
it seems somewhat natural that the theory and formalism of causal Volterra series 
reached its highest development due to the work of (mathematically oriented) 
electrical engineers. This is probably best expressed by the book written by Rugh 
ifTTl . We are now going to apply this previous research to analytic mechanics. We 
will see that enforcing causality in the Lagrangians will result in causal equations 
of motion. By now we should see that the key to deriving purely causal equations 
of motion is to find the correct kernels for our Volterra series. What we will do in 
this section is introduce a new family of kernels to be used in our Volterra series. 
We then examine these kernels to model the nonconservative harmonic oscillator. 
Next, a driving force is added on. Finally, higher ordered potentials are dealt with. 

First, let us define the symmetric and stationary kernels 




(17) 



(V(r - ri)2 + . . . + (r - r„)2 ) 



A-l 




r(A + n-2) 



(18) 



7 



for n >2, and 




+ (r - r„)2 



) , (19) 



r(A + n- 1) 



for n > 1, where, in (fTSt and ( fT^ . < r for i = 1, . . . , n and 



(It is interesting to compare the distributions in dTSb and (fWt with the examined 
in Q.) The restrictions on the Tj should be viewed as part of the definitions of 
the i?^(T) and the r^ir). Requiring that Tj < r is what will give us our causality. 
At first sight the distributions in dTSl and dT^ appear to have nothing to do with 
the fractional derivative operators $q (r) presented in sectional This will prove 
to not be the case presently. 

The easiest case to deal with is the nonconservative harmonic oscillator, which 
we examine now. Consider the Lagrangian 



where < 7 < 2 and C is a constant. Remember that the Lagrangian is a r- 
dependent Volterra series that will be treated as a distribution. From (OTT) we form 
the action 



In we take q(Ti) = for U < a. This implies that C[q; r] = for r < a 
in (OTT) since the Lagrangian in (OTl is causal. Also, the action in (E^ is a t- 
dependent Volterra series. Remember that we think of V[q; t] in (E^ as being the 
anti-derivative of C[q; r] in (OTT) . 

In (OTT) and we should think of r as being 'now'. When we derive our 
equations of motion shortly, we will vary the q{Ti) in (E^ . When we do this, we 
always consider the perturbations of q{Ti) as happening 'now', i.e.. 





(22) 



q{Ti) + hS{T- Ti) 



(23) 
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where h is infinitesimal. Using (E^ in (E^ gives us 



y^l^' J = [ Imr^'^M + mCr^'^M + muj'^r^(T)] -k q dr 

mV Ja 

= I [mj^l [q] + mCJ^l [q] + mu\Jil [q] ] dr. (24) 

We will require that (|24ll is the zero distribution for t > a. This means that 

m^Bl[q]+mC^B:i[q] + mu\B'Aq] = (25) 

for T > a. Equation (l25t is the nonconservative harmonic oscillator's equation of 
motion. The above can also be extended to R2^{t) where A > 2. So, derivatives 
of arbitrary order of g(r) can be included in (|2T]). 

So how did the restriction Tj < r result in the causality in (l25t ? Consider the 
quantity 



S[q;t] = j R^\r)q{ri)6{r - r2) dr2dTidT 

= j Y[-X) (V {t - ^lY + - ^2)^ ) q{Tl)5{T - Ts) dT2dTidT 



(26) 



Since ri < r, we can rewrite (E^ as 



a J a 

t 



r(-A) 

[ $l,(r) * g(r) dr (27) 

where g(r) = for r < a. Hence, the restrictions on the allows us to have the 
^D^ [q] arise from the use of the R2^{t) in the Lagrangians. 

Driving forces need to be handled somewhat differently than the potential or 
kinetic energy terms. This is perhaps not surprising since a driving term is, in 
a sense, outside the 'universe' we are modelling. Let us consider the g(r) as 
our system that moves through some environment described by the kernels in our 
Volterra series. So the Volterra series includes both our system and environment. 
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The driving term is neither of these. We could of course expand our view of what 
the environment or system are to include those mechanisms that give rise to the 
driving force. This, however, may overly complicate matters. Instead, we will 
proceed as follows for, e.g., the driven nonconservative harmonic oscillator. Let 
our new Lagrangians be given by 

C:\q-r] = Cnokr] - r f{n)q{n)dn. (28) 

J a 

Things proceed as before. So our action is given by 

V'[q;t] = VHo[q;t]- [ I f{n)q{n)dndT. (29) 

J a J a 



'f{r)dr. (30) 



Perturbing q(Ti) by h5{T — Ti) results in 

5V'M ^ 5V[q-t] 
5q{T) 6q{T) 

Requiring (l30b to be the zero distribution for t > a gives us 

m^Bl[q] + mC^B^AQ] + mco\B'^[q] = /(r) (31) 

for T > a. So driving forces are easily included in our equations of motion via the 
term 

- r f{n)q{n)dn (32) 

J a 

in our Lagrangians. 

Now we turn our attention to higher ordered potentials. Consider a term like 

Cn[q;r] = ^R>^{T)^q^ (33) 

in our Lagrangians, where n > 2. The term in the action resulting from (l33t is 
given by 

Vn[q;t\ = - f Rl{r)^q^dT. (34) 
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Perturbing q{Ti) as before gives us 



6q{i 



rti{r)^q''~UT (35) 



a Jo Jo 



r / ^ 



r(A + n-2) 



X - n) ■ ■ ■ ^(t- - r„_i) dTn-idTidr (36) 

where (1^ follows from (l35b by the change of variables — > t — Ti. Now we 
switch to spherical coordinates so that ( BHt becomes 



Sq{T) J a Jo r(A + n-2) C„-i 

= / / F7YT 7r,Q{r,r)drdr (37) 

A Jo r(A + n-2) 

where Q(r, r) is that part of the integral in (l36b that depends on r and r after 
integrating over the angles uji,i = 1, . . . ,n — 1. (Note that the restrictions < t 
become restrictions on the cui when we switch to spherical coordinates. Also, the 
exact form of Q{r, r) is unimportant for our purposes, as we now show.) Rewrite 
dTTl as 

^llfyl = fjyU_,{r)Q{r^r) drdr. (38) 

For the potential energy terms we take A = 2 — ri. Then (l38t becomes 

5Vn[q-^t\ 



5q{T) 



Q{0,t) dr 



g"-'(r) dr. (39) 
From (0^ it follows that (B^ . with A = 2 — n, can be rewritten as 

/:n[g;r] = ^ 5 (v/(t - Ti)2 + . . . + (r - t,)^ ) ★ g(ri) ■ ■ ■ g(T„) (40) 
so that (OH) becomes 

V„[g;t] = - f d(Vir-ny + --- + {T-Tny)^ (41) 

lin) ■ ■■qi'Tn) dr. 
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All our potential energy terms in the Lagrangian will be as in (BUt . 

So now we can include driving forces, higher ordered potentials and derivatives 
of arbitrary order in our Lagrangians. We are also able to derive purely causal 
equations of motion by using the kernels in (fTSt in our Lagrangians. A general 
Lagrangian will be given by 




m=2 

where the cj and are constants and < ai < . . . < «„. Allowing aj < 
results in a fractional diffeo-integral equation of motion. 



5 Discussion 

We have presented a rather general method to model nonconservative systems. 
There are two key observations that allowed this to be possible. The main point 
was to find the correct form for the Volterra series expansion of the Lagrangians 
£[g;r]. This form is given in (H^ . In particular, allowing the Volterra series 
kernels to depend on the 'extra' parameter r allowed us to get around Bauer's 
corollary [3| and the result of Dreisigmeyer and Young in Q. Because of this 
expanded definition of the allowed kernels, we were able to find the -R;^(t) that met 
our requirement for a causal equation of motion. That is, we abandoned using the 
fractional derivative kernels $^ (r) in our Lagrangians and actions, as in 1 6 1 . This 
is the second observation that allowed us to achieve our goal of a single retarded 
equation of motion. 

There are a few areas of future research that would be interesting to pursue in 
relation to this paper. Let us initially list the more mathematical of these. First, are 
the -R^(t) the only kernels that allow us to derive the correct equations of motion? 
We have not proved that they are in this paper so, there may be other, more 
appropriate kernels to use. Also, it would be nice to have a rigorous exploration 
of the 'generalized convolution' operation that appears in the C[q] r] and V[q; t]. 
While we treated both the Lagrangians and the actions as distributions, this still 
needs to be given a firm mathematical foundation. Finally, we have only let A G M 
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and n = 2, 3, . . . inthe -R^(r). In fact, forn 7^ 2 only A = was allowed. It would 
be interesting to extend this to arbitrary A G M and n > 2. Allowing A G M for 
any n = 2, 3, ... is relatively straightforward. Allowing for fractional dimensions, 
e-g-, n>2 with n G M, is a more difficult and interesting question. 

For physics, the central problem is probably to develop a Hamiltonian formu- 
lation that corresponds to our Lagrangian one. Also, it would be interesting to 
examine what effects special relativity may have in our formalism. Finally, the 
question of quantization, via, say, the path integral, is also open. All of these areas 
of research may cast light on the question of whether the -R^(t) are the correct, or 
only, kernels to use in our Lagrangians. 
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